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ENERGY ADDITION TO AN ATOMIC HYDROGEN PLASMA 
AT OFF-RESONANT CONDITIONS 
by Richard R. Woollett 
Lewis Research Center 

SUMMARY 

The efficiency of energy addition to a plasma by means of a Stix coil was determined 
for off-resonant conditions. The results included effects of cyclotron damping and axial 
temperature. The coupling resonance was found to be much sharper than had been pre- 
dicted previously. Additional maxima in power absorption were found at slightly higher 
magnetic fields. These off-resonant maxima exhibited magnitudes and breadths that 
made them appear well suited to plasma heating. These new maxima appeared because 
of a beat phenomenon between the forced and natural plasma waves. 


INTRODUCTION 

At present, there is considerable interest (ref. 1) in heating a magnetoplasma by the 
addition of electromagnetic wave energy in the neighborhood of the ion cyclotron fre- 
quency. The calculated magnitude and half breadth of the resonant absorption presented 
by Stix in reference 2 indicate that this technique is an efficient way to heat a plasma. 
However, several problems plaguing the experimentalist make it difficult to verify the 
theory. First, diagnostic techniques for plasmas are often of questionable value. 

Second, in the region of the ion cyclotron frequency, there are theoretically possible 
competing phenomena that could absorb the excited wave energy. An example of such a 
process is the particle resonance associated with atomic ions or with molecular ion-atom 
combinations (ref. 3). Further complications result from the shift in the resonance re- 
sulting from collisions (ref. 4). It would be helpful if theory predicted more than just 
the location of a resonant point. In fact, if the calculated energy addition had distinguish- 
ing features other than a simple resonant maximum, there would be strong evidence of 
verification from just the experimental duplication of the shape of the power curve. The 
present calculation of energy addition by the excitation of the natural plasma waves 



(using a Stix coil, ref. 2) exhibits such distinguishing features. 

The previous analysis (ref. 2) assumed a uniform, fully ionized plasma immersed in 
a large steady-state magnetic field. In addition, it was assumed that there were at most 
only a few electron-ion collisions as the ions traversed under the Stix coil. In the deri- 
vation, however, certain approximations were made in order to present the solutions in 
a compact form. This resulted in expressions only good for a coupling resonant condi- 
tion, that is, where the wave number associated with the radiating coil system was 
identical to the wave number of the natural modes of the plasma. In addition, these 
resonant results were good over limited regions of plasma conditions. Because of these 
restrictions, effects associated with axial temperature and cyclotron damping were lost. 
The theory presented was mainly applicable to a situation where the current boundary at 
the coil had a pure sinusoidal wave form. In this derivation, certain questionable ap- 
proximations were incorporated. This difficulty was later circumvented (ref. 5) by ap- 
proaching the problem in a somewhat different manner. However, axial temperature ef- 
fects, cyclotron damping, and particle energy, which are of interest in the present study, 
are still not covered in reference 5 for a Stix coil. In the present analysis, only a single 
impressed current wave shall be considered. The model used will be the same as the 
one used in reference 2 except that in the present case certain of the previous approxi- 
mations will not be made. The equations will then be applicable at off resonance. This 
procedure will also extend the use of the equations to cover a continuous range of plasma 
conditions that include cyclotron damping and axial temperature effects when significant. 

The results of the analysis are presented in the form of so-called U-factors, which 
represent (as in ref. 2) a nondime ns ional power expression. The nondimensionalizing 
factor is the magnetic field energy associated with a Stix coil for a no-plasma condition. 
The U-factors for both wave energy and the transverse kinetic energy of the particles 
are presented as functions of the steady-state magnetic field for the range from 4000 to 
6000 gauss (the range of interest in current Lewis Research Center experiments). Cal- 
culations are presented for an atomic hydrogen plasma having the densities of 10**, 
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10 , and 10 ions per cubic centimeter and having various axial temperatures between 

the limits of 10^ to 10^ °K. Data for coil current wave lengths of 20, 30, 40, and 50 
centimeters are also presented. A numerical comparison is also made between the 
various dispersion relations involved in the analysis. 


SYMBOLS 


A atomic number of ion 
A ± defined by eq. (4) 
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a m’ b m 
B 


°Q 

B Q (res) 


E 

V z > 

r ± 

J n (z) 

f* 

V z ) 

k 

L 

l 

M 

m 

N 

n 

P 


P± 

Q 

R 


coefficients of Fourier-Bessel expansion for electric field 

nondimensional amplitude of perturbation magnetic field (perturbation field in 
G divided by Bq) 

zero-order axial magnetic field strength, G 

value of magnetic field for which coupling resonance of ref, 2 occurs for a 
given frequency, G 

speed of light, cm/sec 

amplitude of electric field, esu/BQ 

kinetic energy of ions, eq. (10a) 

defined by eq. (32) 

modified Bessel function of first kind and order n (ref. 6) 
defined by eqs. (5a) and (5b) 

Bessel functions of order n (ref. 6) 

amplitude of sheet current, (abs amp/cm)(47r/B 0 ) 

modified Bessel function of second kind and order n (ref. 6) 

wave number, rad /cm 

Boltzmann's constant 

nondimensional coil length, Ico^/c 

coil length, cm 

defined by eq. (25) 

mass of ions, g 

number of wavelengths per coil length (N = 2 for all calculations) 

o 

particle density, ions /cm 

average power transferred from coil to plasma for given sinusoidal coil cur- 
rent, ergs/sec 

peak radio-frequency magnetic field energy in a vacuum (times o>), defined 
by eq. (18) 

(1 ± fi)(/cs 0 ) -1 

storage factor of coil 

nondimensional radius, ro> ■ /c 

Cl 
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<&,£ defined by eqs. (19b) and (19d) 

r radius, cm 

r Q plasma radius, cm 

S(5) defined by eq. (35a) 

Sq, Sj defined by eq. (35a) 

Sq nondimensional axial velocity, (2&T\/mc ^)l/2 

T nondimensional time, to> c ^ 

Tt axial temperature of plasma 

Tq time ion starts path under coil 

t time, sec 

U power absorption factor (U-factor), P/P D 

U Q ,Ui real part of Bessel function of orders 0 and 1, respectively 
VgjVj imaginary part of Bessel function of orders 0 and 1, respectively 
W defined by eq. (12) 

w ci 

Z nondimensional axial distance, z 

charge number of ion 
Zq location at end of coil 

z axial coordinate, cm 


A 

z 

axial unit vector 


a 

477nmc^/Bg 


p 

defined by eq. (8b) 



defined by eq. (31) 


y 

defined by eq. (33) 
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defined by eq. (35b) 


e 

defined by eq. (26a) 
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defined by eq. (26b) 


V 

defined by eqs. (22), 

(23), and (24) 

K 

nondimensional axial 

wave number, k 

a; . 

*D 

(k 2 - 0^1/2 

Cl 
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X 


axial wavelength, cm 
ju. defined by eq. (8b) 

v nondimensional radial wave number 

Vq real part of v 

i/j imaginary part of v 

£ defined by eq. (26b) 

p defined by eq. (8b) 

o ’ 2 axial conductivity, eq. (2) 

cr 2 transverse ‘conductivity, eq. (3) 

t defined by eq. (8b) 

T,$,\l' defined on pp. 19 and 20 
cp defined by eq. (16b) 

nondimensional frequency, co/w . 

O(res) frequency for B Q (res) 

nondimensional frequency oo/oj^, which originates from the vacuum dis- 
placement current term in Maxwell's equation 

a> frequency, radians/sec 

u> c j ion cyclotron frequency, BQe/mc 

Subscripts: 

c cold plasma 

m various natural modes 

p particle 

r, 9, z cylindrical coordinates 

w wave 

x, y, z rectangular coordinates 

1 component of vector perpendicular to the z-axis 

Superscript: 

amplitude factor of the perturbation e i ^ £Z+w ^ 
i designates derivative of Bessel function with respect to its argument 
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THEORY 

Dispersion Relation 

The mathematical model used to develop the equations to determine the power 
transfer is essentially the same as the model used by Stix. Both models used the fol- 
lowing set of assumptions: 

(1) The plasma is fully ionized. 

9 

(2) The I R losses, gravity, viscosity, and electron inertia are neglected. 

(3) The ion temperature is nonisotropic, that is, zero in the transverse direction. 

(4) The distribution of the axial velocity of the ions is Maxwellian. 

(5) All oscillations are periodic in axial distance and in time (i. e. , contain the 
factor e i<kz+tt,t) ). 

(6) The perturbation velocities are small amplitude oscillations perpendicular to a 
steady-state magnetic field. 

(7) Perturbations of the electron and ion axial velocities are ignored. Consequently, 
the derivations will follow very closely those presented in reference 2. 

The essential character of electromagnetic wave propagation through a plasma is 
represented in the dispersion relation. This expression is entirely independent of 
plasma size and depends only on the various transport properties of the gas. Since the 
dispersion relation is completely dependent on plasma currents, which in turn depend on 
particle (ion and electron) velocity, the initial step in the derivation of the equations is to 
obtain relations for the various particle velocities. This is done by first dividing the 
plasma into a massless electron fluid and into constituent streams of ions. Each con- 
stituent stream consists of those ions of the Maxwellian distribution that have the same 
axial velocity. The expression for the average transverse ion velocity is now obtained 
by solving the equation of motion of each constituent stream and adding the various 
streams together. By a similar technique it can be shown that the electron velocity is 
comprised of two parts - one originating from the free motion along the steady-state 
magnetic field and the other from an E X B drift motion. With these expressions for 
the average velocity of the particles under the coil, a plasma current may be derived, 
which when used with Maxwell’s equation results in the following wave equation: 


Jv x (V X E ) - fipE = a jE + icrgE x z] 

(1) 

CT 1 = «(I + +!_-!) + i a (A + - A ) 

(2) 

cr 2 = a(0 + I - I_) + ia(A + + A_) 

(3) 


6 



Equation (1) is the same as equation (7) of reference 2 except for the retention of the dis- 
placement current term Including the displacement currents extends the range of 

validity of the solutions to very high frequencies, such as those corresponding to ion 
cyclotron resonance at high magnetic fields. The quantity £2^ that appears throughout 
the report is identical to £2 but is given the subscript D when it appears in terms that 
originate from the displacement currents. The expressions used for A + and A_ are 
the long-time asymptotic forms presented in reference 2, which are equal to 


XB 0 Z i 

± 9. 47A 1//2 T. 1//2 


e" p ± 


( 4 ) 


The quantities I can be expressed in the form 


K S 


0 




2 /’ P ± 

I, = — e” P± / e x2 dx 


(5a) 


For p ± > 12, equation (5a) is represented to good accuracy by 


i ■ ._S_( 1+ -L * 

* 2p l 2 pI 


(5b) 


Equation (5b) is an asymptotic form of equation (5a) and is obtained by evaluating 


•'0 


y 2 

e x dx for large arguments by L*Hospital rule. For p = 12, the difference 


in the two expressions is less than one in the seventh place. This accuracy is the same 
order as that involved in the evaluation of equation (5a) itself. To evaluate equation (5a), 
a table obtained from reference 7 was used. 

The solution of equation (1) shows that, for a uniform axially symmetric plasma 
where cr^ and Og are independent of position, E r and Eg vary as the unmodified 
Bessel function of the first kind and of the first order. The radial wave number v that 
appears in the Bessel function J 1 (vR) can be written as 




, 2,2 

( CT 1 " k D ) 


a l " K D 


(6) 
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where 


D 



This is the dispersion relation for the propagation of an electromagnetic wave through the 
plasma model described previously. It includes effects of cyclotron damping and vacuum 
displacement currents. 

To determine the real and imaginary parts of v from equation (6), let 

v = Vq + iv^ 

Then 

v 2 = (i'o + i^) 2 = (, 2 - v\) + 2l "o l, l i 


Consequently, 


„ 2 2 2 
Re v = Vq - i^i 

2 

Im = 2vqV^ 

Expressions for Vq and can readily be obtained from these equations as functions 
2 2 

of Re v and Im v , which in turn can be obtained by substituting equations (2) and (3) 
into equation (6) and by collecting the real and imaginary parts of the resulting expres- 
sion. This results in 

R 2 mRm 2 - j) - (T 2 - p 2 )] + 2 t[mt - flp] 

2 2 

jU + T 

Im j,2 = tRt 2 - P 2 ) - (P 2 ~ fl ] + 2 mRxt - /3p] 


( 7 ) 

(8a) 


where 
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! 


fi, = Re 





/3 = Re a 2 
r = Ira (Tj 
p = Im u 2 


(8b) 


With the previous expressions for Vq and it is possible to obtain the radial distribu- 
tion of Eg and E r for a forced axial variation of the electric fields. From these dis- 
tributions it is possible to calculate the power absorbed by the plasma. 


Power Absorption U-Factor 


There are several ways in which the plasma can absorb energy. One of these is by 
means of the kinetic energy of the particles involved, and the other is by a cooperative 
effect of the ions and electrons such that the electric and magnetic field intensities 
within the plasma increase above the vacuum values. The two sets of calculations to 
determine the absorbed energy start from two different and not necessarily self- 
consistent models. 

So far the relations developed do not involve any plasma boundary conditions. When 
energy calculations are made, however, boundary conditions have to be considered. For 
the present calculations, energy is radiated from a Stix coil and is absorbed by the 
plasma. It is then propagated out from under the coil either by a plasma wave or 
through the kinetic energy tied up with the ion motion. The results are presented in 
terms of a relative power efficiency, called the U-factor, and defined by 



(9) 


where P is the average power transferred from the coil to the plasma for a given 
sinusoidal coil current, and is the product of <x> and the peak radio-frequency 
magnetic energy stored in the field of the coil in the absence of a plasma. Actually there 
are other factors (see ref. 2) involved in determining an overall efficiency, but they all 
represent power losses originating from the geometry of the coil system. These addi- 
tional factors are neglected, since the primary concern herein is the ability of a plasma 
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to absorb power, which is solely dependent on the U-f actor. If the plasma fills the 
volume under the coil, the overall efficiency (ratio of power input to the plasma to the 
total power used) is QU/(1 + QU) where Q is the vacuum Q of the coil. 

Power absorption by particle resonance (U ). - Considered first is the energy ab- 

_ r 

sorbed by the particles themselves. It has been shown (ref. 2) that if the transverse 
energy of an ion averaged over an oscillatory cycle is further averaged over a 
Maxwellian distribution of axial velocities, the resulting average energy as a function of 
time (ref. 2) is 


S i ~ me' 


A. 1 


F n + (T-T n ) 


0 ' 2 



2 


A - 

N 

F o + (T - Vy 

|E x -iE y | 2 * 

J 


(10a) 


where Fq is a term independent of time. The power associated with such a phenomenon 
can be obtained by taking the time derivative of this energy: 


P 

P 



dV 


(10b) 


where the integration extends over the volume under the Stix coil. When equation (10a) 
is substituted into equation (10b) 


But 


V 


w .nmc 

Cl 



E + iE 
x y 


2 dV + A 




E v ± iE 


y 




so that 
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( 11 ) 


V 


jrjco C i nmc ^ 


E 

1+ ^ 


iE, 


A + + 


E 

1 - -s£ 
iE, 


¥ 


A - / r |E 0 | dr 


•where 

l = NX 

The solution of equation (1) with the boundary conditions for an infinitely long cylin- 
drical plasma bounded by an azimuthal sheet current having a sinusoidal axial variation, 
that is, j* e *( KZ+wt ), yields 

E e = WJ^yR) 


where 


W = 


j*mR 0 


j ^Rq) k i(' c d r o) 

vR~ : - Kt-sR, 


°J 1 (^R 0 ) D °k 1 (/c d r 0 ) 


( 12 ) 


Jj^Rq) 


This expression can be used with real, imaginary, or complex v. The integral in equa- 
tion (11) now becomes 



(13) 


When v is real, equation (13) yields 



RJj(z,R)dR 


Rn 


r’2 


(J'Rq) 



J^Rq) 


(14) 


When v is imaginary, the integral in equation (13) may be evaluated by setting 
Jj(i^^RQ) = il^jRQ). The result is 
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f 


R E 0 * dR cc 


r °R|j 1 (iy 1 R 0 )| 2 dR= / 0 RI 2 ( yi R 0 )dR 

Jo 


R 0 I f<‘'l R 0> 


2 A 2 

V 1 R 0 


9 I 0 ^1 R 0 ) f 

<‘'i R o> i 

i 2 (kjR 0 ) vw 


( 15 ) 


The remaining case for v complex is somewhat more involved since J 1 (^R) itself is a 
complex number. The integral can be expressed as 

R I J 1 1 2 dr = - 2 ^-j- |[u 1 ( J ,R 0 )U 0 (^Ro) + V 1 (^R 0 )V 0 (t'R 0 )J sec (p 

- [^(i/RqJV^vRq) - U 1 (yR 0 )V 0 (vR 0 )J esc <pj (16a) 


where 


> 


(16b) 


t p = arc tan 


i 2 2 

v\ = (*/g + i/J) 


J 0 (^ R 0 ) = u o(^ R 0 ) + iV o^ R o) 

J^i/Rq) = U^i/Rq) + iV 1 (^'R 0 ) 


The term in equation (11) involving the ratio of electric fields may be evaluated from 
the r component of the wave equation (eq. (1)). Thus, 


E_ 


iEn 2 
d « D - 


( 17 ) 


Consequently, 


12 


E 

1 ±~ 

iE, 


(M t/3) 2 + (t Tp) 2 

2 2 
fi + T 


where p, /3, r, and p are defined by equations (8b). 

Before the U-factor can be obtained, an expression must be derived for the 
quantity P D (eq. (9)). This expression can be given as 

P D =— f B? | B | 2 dV 
D 16ir -Vol 0 


( 18 ) 


where the integration is from zero to infinity radially and is bounded by the ends of the 
coil axially. The magnetic field components under an ideal coil are given by 




B e =0 


= 1<R ~ W*) 


(19a) 


where 


(ft = ij* ORqK 1 (/c d Rq) 

Outside of the ideal coil they are given by 


(19b) 


B„ = g - K.(k d R) 
r S2 1 u 


B& = 0 


\ K 0 < k d r ) 


■\ 


(19c) 


where 
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wv 


(19d) 


$=a ft 

k i(*d r o> 

Adding all the field energies present results in 

P D - 5 

When equations (11) to (14) and (20) are substituted into U p = P p /P D , 

U MO?? 

P ‘ 2K 1 (k d R 0 )I 1 <* d R 0 ) 


where 


V = 


y 0 R 0 


J oW 

'iW 


J 0 (t '0 R 0 ) 2 2 
2p 0 R 0 T + 0 R 0 


J iW 


J 0 W k o("d r o) 

I'nRn + /C~R n 

0 °Ji('’o R o) °Ki((c d r 0 ) 


for v^ = 0 


ri = 


*’i R o 


MW 

I i< 1 'i R o> 


„ V“l R 0* , „ ,2 

- 2u-, R n _ . ~ (i'jRq) 




I i(^i R o) 


"1 R 0 


V»'i R o> 


+ /c d r 0 


K 0 (k d r o) 


for i/g = 0 


?7= R 0 { 


iiW u u k i( k d r o)_ 

[u i (vR 0 )U 0 ( i /R Q ) + V 1 (z/R 0 )V Q (yR 0 )]sec cp- [Uq^RqJV^ Rq) - U^R^V^iy] 

I 


*o + *ff /2 [ u W +v fr' R o>] 


for v complex 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


CSC cp^ 

(24) 
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( 25 ) 


R^U^R^ + £V 1 (t'R Q )" 

k o( k D R o) 

2 

R 0^ U 1^ R 0> - rVi(^R 0 >] 

U^R q ) + V^R 0 ) 

+ “d r o Ki ( Kd r 0 ) 

+ 

U^R q ) + V^R q ) 


£ = VV^O) " U 1 V 0^ R 0) 


k = ^!U 0 (^R 0 ) + ^ 0 V 0^ R 0) 


(26a) 


> (26b) 


J 


With equations (21) to (26a) it is possible to calculate the U-factor due to adding power to 
a plasma by means of the kinetic energy of the ions present. 

Power absorption by means of waves (U w ). - Power absorption by means of waves 

is by far the more important of the two phenomena being considered herein, because the 
magnitude of the U-factor for wave absorption is generally considerably larger except 
perhaps at the coupling resonant condition. This wave energy absorption is accomplished 
by increasing the wave amplitude and changing the phase angle. Before proceeding with 
the analytical formulation of the problem, the boundary conditions must be explicitly 
stated since they play a more integral part in the derivation of wave energy than in parti- 
cle energy. It is considered that the region of the plasma which responds to the forced 
radio frequency wave lies inside a s, box” beneath the Stix coil. This plasma box is 
bounded on the left and right, outside of the Stix coil, by a cylindrical plasma that ex- 
tends to infinity and has a radius r^ equal to that of the coil. The plasma box is also 
bounded by the coil itself. The solutions for the plasma box and the infinite plasma out- 
side the coil are joined by making the wave amplitudes continuous at the boundary. Out- 
side of the coil, the wave motion is assumed to be a natural mode of oscillation of a 
cylindrically bounded plasma in which the wave amplitude is constant. In the left hand 
region, the amplitude of a right running wave is zero. The right running wave squirts 
out the right end of the coil and propagates axially to infinity (ref. 2). An analogous 
phenomenon occurs with the left running wave. Underneath the exciting coil, it has been 
assumed that the forced plasma wave is equal to that generated by a Stix coil of infinite 
length. The forced wave, that is, J^(yR), in the plasma box under the Stix coil can be 
Fourier decomposed into an infinite series of J^z^R); that is, 
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- E a m J l(V R ) 


XL 

where i; is the radial wave number for the m 111 radial mode of oscillation of the free 
cylindrical plasma. The forced electric field is given by 

^ e i(«Z+S2T) _ WJi( „ R)e i(«Z + fiT) 


Consequently, 

00 

’ X Wa m J l ( V R ) 

m=l 


The natural mode relation used was obtained from expressions previously developed 
(ref. 8) for cold collisionless plasmas (for which no cyclotron damping occurs). For a 
given fl, Bq, and n these natural modes can be obtained from the simultaneous solution 
of the following two equations (ref. 8): 


and 



(27) 


v c 


Rr 


J p(*c R 0 ) 

R o) 




(28) 


Equation (27) is the dispersion relation for a cylindrically symmetric plasma where the 

i(/c c Z+S2T) 

perturbations have the form J ^(^ c R)e . Equation (28) is obtained from the wave 

equation and the * boundary condition” for a vacuum boundary of a plasma located at 
R = Rq. The azimuthal electric field, along with its associated magnetic field, were 
both assumed continuous across the boundary. In the derivation of equations (27) and 
(28), the axial electric field was set equal to zero (negligible electron mass). It is not a 
particularly straightforward calculation to find a solution which is designated and 
K m from these expressions. The technique used herein is outlined in appendix A. 

Once the values of the wave numbers associated with the natural modes of plasma 
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Figure 1. - Schematic of spatial variation of wave 
amplitude and coil current 


oscillation have been determined, an expression for 
the electric fields under the coil can be derived. In 
order to satisfy the boundary conditions existing on 
each end of the coil, a natural mode wave must be 
added to the forced disturbance. This is analogous 
to adding a complementary solution of a differential 
equation to its particular solution. The electric field 
under the coil can be represented by 

E e e 1 ( KZ+nT ) = ^ f m R)e 1 <' cZ+nT > 

+ E b m J l<V R ) ei( ' t ” Z+nT) <29> 


where the b m f s are coefficients to be determined by the boundary conditions. (Only a 
term equivalent to the first series on the right hand side of eq. (29) was required for the 
fields in the particle model. ) The boundary condition requires that the electric field for 
each left running radial mode vanishes at the right end of the coil, that is, at Z = Zq 
( see fig. 1). This results in 


Wa e 
m 


i/cZ n i/c,Z n 

°+b e m °=0 
m 


and 


b = -Wa e 
m m 


i(K-V) Z 0 


Consequently, equation (29) become s 


g.^.i(/cZ+ftT) _\ ^m J l^m R ) F _ e " 1 ( /c V^ Z-Z (/| e i 

Z — i 1/2 - v m ^ J 


i(/cZ+S2T) 


(30) 


If it is assumed that the terms of the Fourier Bessel expansion of J^iTt) are all ortho- 
gonal, then the following expression for /3 m can be developed by normal techniques: 
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where 


and 


2i£2j*y 

m Rot 1 + <9 2 - »(-m R 0 )" 2 ] J l<>'m R o) 



J 0^V / 2 0 2\ 


Y = 


Kpfe - n p) 1/2R o) 

- n p) 1/2R o) 


VR r 


Jq(^Rq) 
J ^(i/Rq) 


+ KpRg 


k 0 ( k d r o) 

k i(«d r o) 


( 31 ) 


(32) 


(33) 


This expression for |3 varies somewhat from that of reference 2 because of the in- 
clusion of several terms that were dropped in that reference. The assumption that all 
the J^(t' m R) , s are orthogonal, of course, is incorrect. However, calculations have 
been made of the greatest possible error evolving from this approximation. The results 
indicate that the error in the various coefficients of the expansion is less than 1 percent. 
The expression for y given by equation (33) can still be used when v is complex. It 
can be noticed that when v is either real or imaginary, y is real; when v is complex, 
y is complex. Equations (30) to (33) represent the electric field of a particular plasma 
mode set up by distributed currents within the plasma. If this field is evaluated at the 
coil radius, total power input into the plasma may be obtained by integrating over the 
surface of the plasma boundary the product of the electric field and the current required 
to create it. When this operation is performed, there is an expression for the average 

■fVi 

power going into the m in radial mode (Pj m : 

W HI 
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where 


8(8) = 1 - COS »j +K«>- »»n 28) = Sq + ^ 

2o 2 


(35a) 


and 


6 = (k - « m )L/2 (35b) 

It may be noticed that equations (34) and (35a) are identical to equations (22) and (23) 
of reference 2. The region of O in which damping is important will be shown to be ex- 
tremely small and centered about the point 0=1. If there is no damping present, only 
the real part of S(5) enters in the expression for real power; the imaginary term con- 
sequently represents the reactive power. The previous expression for power has the 
units of ergs per second. If the resistive power (eq. (34)) is nondimensionalized with 
respect to the same expression (eq. (20)) used with the particle power, then the U-factor 
for the resistive wave, m mode, is 


4L(/c - K m )T 


(U ) = 

' w'm 


2r; 


( 2 2 2 V „ 2 2 

("0 - "1 - v m) + 4 Vl 

^~2 __ ^ 


1 + 


L 




(36) 


For Vq or v j = 0, 

T = [ s o(V ,; r ,, i) + 2 Yi s i]’' 

For v complex, 

T - (Re *)[s 0 <* 2 - v\ - + 2( V ,Sj)] - (Im t)^* 2 - v\ - v 2 m ) - 2 V ,S 0 ] 


where 


19 



Ill lllllllllllllllllllllllllllllllll 


Re 


Re $ + f(/c m )j j*Re $ + f(/c)J + (Im #) 2 
[Re $ + f ( K )] 2 + (Im $) 2 


Im ^ = 


Im $[f ( k) - 

[Re $ +f(/c)] 2 + (Im $) 5 


Re $ = 


_ ftoltWo - "l V P> + V 1<^1 U 0 + W] 


vl + v\ 


Im <& = 


_ go^WJp * VV - v^o - •’lV] 


U 2 + V 2 




and where the U’s and V's are defined by equation (16b). 

All the relations required to determine the efficiency of adding energy to a plasma 
by means of the wave energy and the particle energy have now been presented. The ex- 
pressions do not require that calculations be made at a natural wave resonance condition 
Calculations were performed at various ion densities, frequencies, axial ion tempera- 
tures, and plasma radii. They were performed for a fully ionized atomic hydrogen gas. 


RESULTS 

The results will be discussed in four sections, each of which is fairly independent of 
the other. The first section, Wave Number, examines the natural modes of a cylindrical 
plasma for a range of plasma radii, wave frequencies, and particle densities. The 
second section, Absorption Characteristics, discusses, without reference to the actual 
magnitudes of the parameters involved, the trends of energy addition to the plasma. The 
third and fourth sections, "Wave absorption (U-factor)" and "Particle absorption 
(U-factor), " present the numerical results of the energy addition calculations. 
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Wave Number 


Since the natural modes play an integral part in the development of the theory, it is 
of interest to examine their characteristics in detail before the results of the U-factor 
calculations are presented. The natural modes used herein were obtained from the solu- 
tion of the dispersion relation (eq. (27)) and the equation obtained from the boundary con- 
dition (eq. (28)). The analysis technique used is presented in appendix A. A graphic 
presentation of the solution is given in figure 2. Moreover, various trends can be pre- 
dicted from these figures. The dispersion relation and boundary condition solutions are 
presented as solid and dotted curves, respectively. The natural mode, obtained from 
the simultaneous solution of equations (27) and (28) are consequently represented by the co- 
coordinates of the intercepts of the two curves. In figure 2(a) the boundary condition is 
shown for the first three modes of a plasma with a 5-centimeter radius. In figure 2(b) 
the boundary condition is presented for the first mode of plasmas of various radii. 
Qualitatively, a higher mode can be compared to the first mode for a plasma of smaller 
radius. It can be noticed from these figures that the axial wave number k is practi- 
cally independent of the boundary condition while the radial wave number v is practi- 
cally independent of the dispersion relation. In other words, the plasma radius and mode 
determine the radial wave number while the frequency determines the axial wave number. 
The magnetic field value of 4200 gauss in these figures is typical of the values used in 
the Lewis plasma heating experiments (ref. 1). 

A summary of the fundamental mode wave numbers presented as a function of radius 
is shown in figures 3(a) and (b). The insensitivity of /c to radius is again demon- 
strated, as is the dependence of the radial wave number v . Effects of plasma density 
can be seen by comparing the two figures (note that the scale for k is not the same on 
these two figures). 

An upper limit to the radial wave number v m can be found. If the boundary were an 
infinitely conducting surface rather than a vacuum, the radial modes would be deter- 
mined by the zeros of J^(yR). The resulting wave numbers would be larger than the 
ones associated with the vacuum boundary. Moreover, the actual value of the first zero 
is less than This upper limit to the first radial mode is included in figure 3 for 

comparison. 

The dispersion relations presented in figures 2 and 3 however, do not include any 

cyclotron damping effects. An indication of when these effects are significant can be ob- 

2 2 2 2 
tained by comparing the real part of v , that is, - v^, from equation (6) with v 

from equation (27). The dispersion relation of the cyclotron damped wave is plotted in 
figure 4. By comparing figure 2 with figure 4, it can be seen that although the equations 
(27) and (6) start from different physical models (the latter admits a Maxwellian distri- 
bution of axial velocity and the former does not), the calculated results are very nearly 
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Nondimensional 

frequency, 




JNondimen- 
sional fre- 










Axial wave number, k ? 


Figure 4. - Dispersion relation for a forced axial distribution including effect 
of cyclotron damping. Particle density, 10* 1 ions per cubic centimeter; 
magnetic field, 4. 2xl0 3 gauss. 


the same. It can also be noticed that although the quantity approximates the 

nondamped expression closely, there is a difference in the expressions for values of S2 
close to 1. The effect of the damping process becomes more and more dominant as the 


frequency approaches closer and closer to unity. The first major difference for 

T, = 10 4 °K is noted at S2 = 0. 993, and for T. = 10 5 °K at fi = 0. 985. Above these 

1 1 2 2 
frequencies but for £2 still less than 1, the peak of - v ^ rapidly decreases and be- 
comes negative. 


Absorption Characteristics 

Two different power absorption factors were determined: those associated with the 
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Figure 5. - Schematic of general characteristics of wave power 
absorption factor. 


wave energy of the various plasma modes 

(U„ ) and those associated with particle 
' w'm 

energy U . The total wave energy escap- 
ing from the coil is the sum of the energies 
of all the individual modes present. For 
this report, however, only wave energies 
of the first three modes were added. This 
sum is represented as U , so that 

U w = (l Vm=l + (I Vm=2 + ( U w>m=3 


The general characteristics of U w observed from the present calculations are shown as 
a function of Bq in figure 5. The first three peaks measuring from the left represent 
the coupling resonance for each of the three modes. The scale of Bq in this region is 
greatly enlarged (about 50 times). These maxima occur (see eq. (36)) when v Q « v . 
The fourth and fifth peaks in figure 5 do not occur near any resonant conditions and 
therefore are referred to as off-resonant maxima. They originate from that factor of 


the expression (34) represented by 



The peaking is actually a beat phenom- 


enon (ref. 5, p. 98) developed between the forced and natural modes (see eq. (29)) that 
exist under the coil. 

This description differs from that presented in references 2 and 5, wherein it was 
assumed that the resistive plasma loading does not vary appreciably near the wave reso- 
nant condition. When such an assumption is made the general shape of the U w curve is 
similar to the dotted line in figure 5. Consequently, references 2 and 5 predict a wave 
resonant maximum of large breadth whereas the analysis presented herein predicts 
several extremely narrow wave resonances along with a medium-width off-resonant 
maximum. The narrow wave resonant condition is slightly shifted to the right from the 
resonant condition of reference 2 by a Doppler effect. The off-resonant peak, however, 
is considerably displaced from the reference 2 maximum. 

It is not obvious from the equation for the particle Up-factor (eq. (21)) that maxima 
will occur at all Bq fields where some natural- mode axial wavelengths equal the wave- 
length of the coil. However, plotting Up at high T^ does show this characteristic, as 
depicted in figure 6. The three peaks are associated with the three natural wave modes 
indicated in figure 5. The displacement of the maximum from S2 = 1 includes the 
Doppler effect associated with moving particles. Increasing the values of T^ increases 
this Doppler effect and consequently increases the displacement. When S2 » 1, however, 
cyclotron damping may become so strong that the energy associated with particle motion 
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is essentially "damped” out. The term 
cyclotron damping, as in reference 2, is 
associated with the radial flow of energy. 
The damping modifies the plasma currents 
in such a way as to reduce the magnitude 
of the transverse electric fields. Since 
particle velocity depends on the size of 
the fields, kinetic energy will also de- 
crease along with the wave amplitude. At 
these conditions, figure 6 would not be 
typical of the Up plots. Near the couple 
resonant condition the curves of Up are quite similar to the curves of U w ; each gen- 
erally exhibits maxima at the couple resonant points. However at larger magnetic fields 
Up is very small, whereas U w is generally significant. This small value of Up in- 
dicates that very little energy is transferred to the particles at the off-resonant maxima. 

To investigate the U-factors in more detail, the results of the actual calculations 
must be examined. Since the steady-state magnetic field strength is the primary vari- 
able in most cyclotron-resonance plasma experiments, the U-factor plots were calcu- 
lated and are presented as functions of magnetic field. A general description of the flow 
diagram for the program used to calculate the U-factor is presented in appendix B. 

Because of the interest in comparing the present results with previous work, all the 
plots start at a magnetic field consistent with the cold plasma coupling resonant condi- 
tion determined from the following dispersion relation (eq. (26) of ref. 2): 

S2 2 = — — (37) 

a 


2 _ 10 o 

where a//c = o. 5x10" 1 nA , as given by equation (10) of reference 2, and where 
£2 = u>/(BQe/mc). At the coupling resonant condition A, the wavelength of the natural 
mode is equal to the wavelength of the coil. At this condition, the magnetic field will be 

q 

referred to as Bp(res), which is generally set at 4. 2x10 gauss. The frequency given 
by equation (37) for Bp(res) will be referred to as fl(res). Once the magnetic field for 
coupling resonance is selected, the dimensional frequency can be calculated from equa- 
(37). Consequently, the parameters n, BQ(res), l, and f2(res) appearing in each plot 
satisfy equation (37). 

It was noticed during the calculations that quite frequently negative U w ~factors oc- 
curred at magnetic fields somewhat smaller than the fields associated with the couple 
resonant condition. This would suggest a plasma wave radiating energy to the coil and 



Figure 6. - Schematic of general characteristics of particle power 
absorption factor. 
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therefore would require an additional mechanism to the ones being investigated just to 
sustain its existence. This added complexity was not considered herein, and regions of 
negative U-factors were not plotted. 

Wave absorption (U w ~factor). - The U w -factor plots for a coil length of 80 centi- 
- 11 12 14 

meters and for particle densities of 10 ,10 , and 10 ions per cubic centimeter are 

presented in figure 7. Since the phenomena of coupling resonance occur in regions so 
small that meaningful data cannot be obtained from plots including both resonant and off- 
resonant peaks, the data for the two effects are presented separately. Figures 7(a), 

(b-1), and (c-1) present the data for the off-resonant part of the U -factor for particle 

111214 w 

densities of 10 ,10 and 10 , respectively. Figures 7(b-2) and (c-2) present the 

12 14 

data for the coupling resonance for densities of 10 and 10 , respectively. There is 

no resonant curve for 10 ^ since the effect is damped out at these densities and fre- 
quencies. By comparing the location of the maxima of the various curves presented, the 
trends of the peaks as a function of density may be obtained. However, since the S2 
corresponding to resonance changes with the density (eq. (37)) and since BgCres) re- 
mains unchanged in the various curves of figure 7, u> is different for plots with different 
densities. 

The effect of density on the location and the magnitude of the various maxima is 
portrayed more directly in figure 8. Here the dimensional frequency cu is kept con- 
stant and BQ(res) is permitted to vary. In figure 8(a) the location of the first and 
second off-resonant peaks is shown together with the location of the resonant peak as 
given by reference 2. The resonant maximum calculated herein would fall in the narrow 
region between the curve from reference 2 and that for the first off-resonant maximum. 
All three curves are of similar shape, increasing monotonically with density in the 
range considered. It would be difficult to determine which curve was being observed ex- 
perimentally. The difficulty is further aggravated because the half breadth of the res- 
onant peaks is generally much smaller than that of the off-resonant peaks. Conceivably 
the former would not even be observed. Although the half breadth generally cannot be 
determined from the results, the actual width of the resonance is what is significant. In 
almost aH cases investigated (which covered a range of particle densities from 10 ^ to 
10 1 ions/cm 0 , a range of magnetic fields from 10° to 10 G, and various coil lengths of 
40 to 100 cm), the width of the first maximum was of the order of 5 gauss or less. Non- 
uniformities in density and the steady-state magnetic field, along with various dissipa- 
tive phenomena, however, could decrease the magnitude and broaden the width of the 
resonance. With such conditions present, the couple resonance could become observ- 
able. The relative magnitudes of the peaks are shown on figure 8(b). The first and 
second off-resonant peaks exceed in magnitude the resonant maximum from reference 2. 
The resonant maximum calculated herein is not shown even though its magnitude would 
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Wave power absorption, 



Axial magnetic field, Bq, G 


(a) Particle density, 10 11 ions per cubic centi- 
meter; frequency for value of magnetic field 
for which coupling resonance of reference 2 
occurs, 0.99549; off- resonant maxima. 



(b-1) Off-resonant maxima. (b-2) Resonant maxima. 


(b) Particle density, 10 12 ions per cubic centimeter; frequency for value of magnetic field for which coupling resonance of reference 2 
occurs, 0.96037. 



Axial magnetic field, Bq, G 

(c-1) Off-resonant maxima. Resonant maxima. 

(c) Particle density, lO 1 ^ ions per cubic centimeter; frequency for value of magnetic field for which coupling resonance of reference 2 
occurs, 0.33533. 

Figure 7. - Wave U-factor as function of magnetic field for various plasma densities. Magnetic field for which coupling resonance (rf 
reference 2 occurs, 4. 2xl0 3 gauss; coil length, 80 centimeters; plasma radius, 5.0 centimeters; axial temperature of plasma, 10 4 K. 


30 



Wave power absorption, Axial magnetic field, Bq, 



Particle density, n, ions/cm 3 
(b) Magnitude of maxima. 

Figure 8. - Characteristics of maximum efficiency as functions of particle density. Frequency, 3. 86x10^ radians per 
second; coil length, 80 centimeters; plasma radius, 5 centimeters; axial temperature of plasma, 10 4 °K. 
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Wave power absorption, U w Axial magnetic field, Bq, 




(b) Magnitude of maxima. 

Figure 9. - Maximum U-factors as function of coil length. Frequency, 
3. 8x10^ radians per second; particle density^ lO 1 ’ ions per cubic 
centimeter; axial temperature of plasma, ltr 0 K; number of wave- 
lengths per coil length, 2. 


exceed the others. Its previously dis- 
cussed narrowness limits its signifi- 
cance in a real plasma. 

The U-factors for coil lengths of 
40, 60 and 100 centimeters were also 
investigated. Since the general shape 
of the curves is similar to the data 
for a coil length of 80 centimeters, 
however, only the location and magni- 
tude of the off-resonant maximum and 
the couple resonance peaks of refer- 
ence 2 are presented in a summary 
plot (fig. 9). Here again the dimen- 
sional frequency u> is constant at the 
same value as in figure 8. It can be 
seen from figure 9(a) that there is a 
very small shift in the first off- 
resonant maximum as the coil length 
increases. This shift is quite similar 
to that calculated from reference 2. 

A larger effect, though still rela- 
tively small, is the decreasing magni- 
tude of the U w of the first off- 
resonant maximum as the coil length 
increases (fig. 9(b)); the U w ~factor 
drops roughly 30 percent as 7 in- 
creases from 40 to 100 centimeters. 
Although the magnitude of U w calcu- 
lated from reference 2 is considerably 
less than that calculated for the first 
off-resonant maximum, the two 
U w -factors plotted as a function of l 
do parallel each other for the l 
range investigated (fig. 9(b)). 

One of the parameters of the sys- 
tem that is difficult to determine ex- 
perimentally is the plasma radius. 

The radius used to calculate U w up 
to now has been 5 centimeters. To 
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(a) Axial temperature of plasma, 10^ °K. 



(b) Axial temperature of plasma, 1. 2x10^ °K. 

Figure 10. - Wave U-factor of couple resonance as function of magnetic field and axial ion tempera- 
temperature. Particle density, 1012 ions per cubic centimeter; magnetic field for which coupling 
resonance of reference 2 occurs, 4. 2xhP gauss; frequency for value of magnetic fields for which 
coupling resonance of reference 2 occurs, 0.96037; coil length, 80 centimeters; plasma radius, 

5 centimeters. 

ascertain the effect on efficiency, U calculations were also made for plasma radii of 

w q 12 

0. 5, 1. 0, 2. 5, and 10. 0 centimeters at B = 4. 2x10 gauss and n = 10 ions per cubic 
centimeter. The results indicate that as the radius increases from 0. 5 to 10 centi- 
meters, the U w -factor of the first off-resonant maximum increases from 1. 28 to 1. 47. 
Although the effect of a reasonable change in radius is not large, the results do suggest 
that the plasma should be as big as possible for the most efficient power transfer. The 
calculations, of course, assumed that the coil radius remained equal to the plasma 
radius. If the coil radius increases for a fixed plasma radius, there is an additional loss 
in efficiency that has not been considered, that is, the volume effect of tying up field 
energy between the plasma and the coil. This does not contribute in any way to the 
energy added to a plasma but only requires more current in the coil circuit in order to 

compare cases of equal electric fields at the plasma surface. Consequently, there re- 

2 

suits an increased I R loss. 

It can be seen by examining a large number of calculations that the U w -factors in 
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the region of the off-resonant maxima are not 
significantly affected by varying the axial tempera- 
ture; in fact, the magnitude and location of both 
off-resonant maxima remain constant. There is a 
temperature effect, however, on the couple res- 
onance. If such effects are not damped out com- 

Figure li. - Schematic demonstrating Doppler effect pletely, then a variation in streaming velocity of 
on power absorption factor near coupling resonance. * 

the ions will modify the shape and location of the 
couple resonant maximum. This is demonstrated by comparing the first mode maximum 

O 

in figure 10(a), which occurs at roughly 4. 2x10 gauss, with the first mode maximum in 
figure 10(b), which occurs at 4. 22x10 gauss. This shift is similar to the Doppler shift 
depicted in figure 11. The two symmetric peaks for T^ > 0 represent the Doppler fre- 
quency shift corresponding to a> ± kV z ( see eqs. (3) and (4) of ref. 2). Because of the 
strong cyclotron damping effects in the neighborhood of S2 = 1, however, the left mem- 
ber does not always appear. At a density of 10^ where S3 « 0. 996, the damping is so 
strong at Tj = 10 ^ °K that even the right-hand member does not appear. Figure 12 
presents the calculated shift of the natural mode resonance along with a simplified 
Doppler shift which was determined from the expression 

kSginc^ 

en(res) 

It can be seen that the two calculations do not agree. The displacement of the first mode 
coupling resonance was much smaller than that calculated by the simplified expression. 

At 10 5 °K, the displacement was less than one-fifth of that from the simplified equation. 

In any case, an interesting observation may be made. If a cold plasma is to be 
heated to high temperatures at resonant coupling conditions, it may become necessary to 
adjust the system to stay in resonance as the plasma is heated. As pointed out pre- 
viously, operation at an off-resonant peak would be less affected by changing tempera- 
tures. 

Particle absorption (Up -factor). - When the particle absorption factor Up and the 

wave absorption factor U w are compared, it is recalled that the fields used for the two 
factors are determined from different models. 

The variation of the particle energy with respect to the magnetic field is quite simi- 
lar, near the couple resonant condition, to the variation of the wave energy. Whenever a 
coupling resonant maximum appears in U w there is a corresponding maximum in Up. 
Consequently, Up is generally very sensitive to variations in Bq in the vicinity of the 
couple resonant condition. Whenever the coupling resonant effects are negligible or non- 
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K 



10 3 10 4 10 5 10 6 


Axial temperature of plasma, Tj, °K 

Figure 12. - Location of couple resonance as function of axial temperature. Particle density, 10 12 ions per cubic 
centimeter; magnetic field for which coupling resonance of reference 2 occurs, 4. 2x10^ gauss; frequency for 
value of magnetic field for which coupling resonance of reference 2 occurs, 0.96037; coil length, 80 centimeters; 
plasma radius, 5 centimeters. 

existent in one energy, there is generally a corresponding result with the other energy. 
Consequently, whenever the couple resonance wave is totally damped, the particle ab- 
sorption also becomes negligible. 

The particle resonance also depends strongly on the particle axial temperature, just 
as the U w -factor does. These latter results are depicted in figure 13 in a very expanded 
scale for the conditions of BQ(res) = 4. 2x10 ^ gauss^ n = 10^ ions per cubic centimeter, 
and for Tj = 4x10^, 6x10^, 8x10^, 10^, and 1. 2x10 °K. This figure shows that as the 
axial velocity of the ions decreases (lower streaming temperature), the magnitude of the 
maxima increases, the half breadth decreases, and their location moves to smaller 
values of Bq. The half breadths associated with these maxima, however, are generally 
so small that it would be unlikely that the resonance conditions could be experimentally 
detected. It should be reemphasized that any conclusions concerning the half breadth of 
the couple resonance are questionable. Since variations in Bq are inevitable in an ex- 
perimental system (even exceeding the half breadth of the peak), the effectiveness of the 
Up -factor should perhaps be related to the area under the peaks rather than to their am- 
plitude and half breadths. Consequently, such resonances could possibly be detected. 

On the other hand, it would probably be unlikely that the large peak magnitudes could 
actually be realized. 
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for value of magnetic field for which coupling resonance of reference 2 occurs, 0. 96037; coil length, 80 centimeters; 
plasma radius, 5 centimeters. 

At magnetic fields larger than the couple resonant condition, including the fields for 

the off-resonant peaks of U w , the particle absorption factor U p is generally so small 

in comparison as to be negligible. One exception to this generalization was found in the 

calculations. This case, which is shown in figure 14, indicates the variation of the 

1 2 

U-factors as a function of Bq for a 40-centimeter coil (X = 20 cm) at n = 10 ions per 
cubic centimeter, BQ(res) = 4. 2x10** gauss, r = 5 centimeters, and = 10^ °K. Under 
these conditions, the location of the first off-resonant peaks for U w nearly coincides 
with the location of the couple resonance. This coincidence results principally from the 
reduction in coil length. As was shown in figure 9, the off-resonant peak moves to lower 
Bq fields as the coil length decreases. On the other hand, the coupling resonance for 
waves at a given temperature shifts to higher B fields as may be seen from the simpli- 
fied Doppler shift equation. For the case shown in figure 14, both peaks have shifted the 
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Figure 14. - U-Factor as function of 
magnetic field. Particle density, 

10 12 ions per cubic centimeter; mag- 
netic field for which coupling reso- 
nance of reference 2 occurs, 4. 2x10^ 
gauss; coil length, 40 centimeters; 
plasma radius, 5 centimeters; axial 
temperature of plasma, 10 4 °K. 


proper amount so that they will coincide. It should be noted that both absorption peaks 
are relatively narrow (compare with fig. 7(b)). However, the possibility exists that for 
some conditions a significant amount of energy could be added to the plasma directly as 
particle motion, without going through the two-step process of generating and then 
thermalizing a wave. 


CONCLUSIONS 

The model used herein to calculate the efficiency of heating a magnetoplasma by the 
use of a Stix coil is essentially the same as that used in reference 2. However, certain 
approximations of reference 2 restrict the use of the equation to the coupling resonant 
condition. In the present case, these approximations are not made. Consequently, the 
relations developed are applicable off resonance. 

The actual power transfer is accomplished by two concomitant processes, the most 
important being wave generation. Here the Stix coil radiates energy to a cylindrical 
plasma core, which in turn transmits the wave energy out from under the coil in an axial 
direction. Accompanying the wave process is an energy addition to the streaming ions. 
The transverse particle energy is necessary in order to establish the traveling axial 
waves. When these powers are nondimensionalized and are plotted as functions of the 
superimposed magnetic field, the resulting curves have sufficient complexity to permit a 
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more significant comparison with experiment than just the location of a resonance. The 
following seven characteristics were common to all the response curves. 

1. When the power absorption factor (i. e. , U-factor) for waves is plotted as a func- 
tion of the magnetic field, there always exists an off-resonant maximum that appears 
close to the couple resonant point of the first mode excitation. These additional maxima 
originate from a beat phenomenon between the forced and natural plasma modes. This 
new maximum is located on the high magnetic field side of the couple resonant condition. 
Since this off-resonant peaking has a considerably larger half breadth than the couple 
resonant half breadth and since it lies so close to it, a possibility exists that it may be 
mistaken for the coupling condition. The location of the off-resonant peak behaves in a 
manner similar to the natural mode resonant condition; consequently, a mistaken identity 
can easily be made which could lead to estimates of density that are large. This descrip- 
tion differs from that presented in references 2 and 5 since in these references it was 
assumed that the resistive plasma loading does not vary appreciably near the wave res- 
onant condition. When such an assumption is made, the wave U-factor monotonically 
decreases in the region between the wave resonant and the off-resonant maximums, 
which results in a prediction of a single large breadth maximum. 

2. An additional maximum also appeared at magnetic fields larger than the first off- 
resonant maximum fields. It also originates from a beat phenomenon. The half breadth 
for this peak, however, is considerably larger than the first off-resonant maximum. Its 
amplitude is roughly one-third that of the former. 

3. Both the first and second off-resonant peaks are insensitive to variations in the 
plasma size and axial temperature. However, the first off-resonant U-factors for 
waves did show an inverse wave length sensitivity for the shorter wavelengths. 

4. The magnitude of the first and second off-resonant peaks was within an order of 
magnitude of the couple resonant maxima of reference 2; the first off -resonant peak was 
always larger. 

5. The half breadth of the natural wave resonant maximum was found to be small, so 
small that it could be difficult to experimentally detect for coil wavelengths greater than 
30 centimeters. However, because of nonuniformities in the experimental densities and 
fields, it is possible that such peaks could be broader than indicated herein. At a small 
density, where such maxima would be very near the ion cyclotron condition, the couple 
resonance is completely damped out. As density increases and the resonance moves 
away from the ion cyclotron condition, less cyclotron damping exists and the maximum 
is observed. 

6. The maxima of particle energy addition occur at the same magnetic field as those 
of wave energy addition for the natural wave resonances. Particle energy was generally 
negligible at all other values of field strength even at values where the wave energy 
peaked at the off-resonant maxima. For a short coil case (coil length, 40 cm), however, 
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the particle energy was found to be the same order of magnitude as the wave energy at 
the first off-resonant maximum. This comes about because the couple resonant and the 
first off -resonant maximum occur at the same field strength. 

7. As the axial temperature increases, the ’’coupling” maxima of the particle 
U-factor move to larger fields. The half breadth increases and the amplitude of the 
maximum decreases during this variation. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 2, 1965. 
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APPENDIX A 


SOLUTION OF EQUATIONS (27) AND (28) 

Ordinarily a Newton-Raphson technique would be used to solve equations (27) and 
(28) simultaneously. However, because the solution of the equations falls extremely 
close to the singularity of equation (27), that is, 

jo* + (k c - Op)J + ( K c ~ ^p) = ® 


the correction interval for the first guess could place the second guess on another branch 
of the curve of equation (27), which would result in a hopeless condition. This difficulty, 
however, could be circumvented by using a technique that may be considered as an exten- 
sion of the normal Newton-Raphson equations. 

If two new functions i// and cp are defined using equations (27) and (28) 


> K r ) = 


-fiV + fi 2 [2a + (k 2 c - flp) ](*c - n p ) “ ( K c -°D ) 2 2 

+ (« 2 * ° d )] + ( K l ' 


^c 


2 n ,^ H gofej - $ „ _ J fl». R 0 > 

c w "o , , \l/2 \ c Jj(» c Ho) 


K l(( K c - 4) R o) 


(Al) 


and if these are expanded in Taylor series 


Si// Si// 1 2 9 2 ^ . . 9 2 i// 1 A 2 9 2 ^ 

il/(v + Av, K + Ak) = \p + Av —+ An — +- Av — ~ + Av Ak -t- - A/< — — 


dv Sac 2 


Si/ 


Si/ 3/c 2 g ^2 


(A2a) 


/ . x . 89 ? A dcp 1 A 2 9 2 ^ » . <P 1 . 2 9 2 <p 

<p(i/ + Ai/, K + A k) = cp + Av c + A/c — + — Av — r + Ai> An + — A/C 

dv Sac 2 g ^2 dvdn 2 gf£ 2 


(A2b) 
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Jn(x) 


/rRegions of 
/ \ possible solutions-^ 

/i /» 

/ t / ‘ 

/ \ / \ 


A solution exists when 



Figure 15. - Schematic of zero- and first- 
order Bessel functions. 


\J/(v + Av, k + A k) = 0 


cp(v + Av, k + A k) = 0 


When this latter condition is used, equations (A2) re- 
duce to two algebraic quadratic equations with two un- 
knowns Av and An. Eliminating one of the unkonwns, say Av, results in a quartic 
equation in A k which must be solved. This can be done quite readily with the use of any 
one of a number of machine programs. Of the four roots that exist, the one with the 
smallest absolute value was used in the determination of the correction interval. The 
values of v Q and k c for a solution are designated v m and « m . The only remaining 
problem is that of obtaining the first guess of the solution. Equation (28) is some help; 
since the right-hand side of equation (28) is always negative, the sign of mus ^ 

always be opposite to that of Ji(y m R n ). Because of the sinusoidal nature of the Bessel 
functions, the regions of possible solution are values of x between the n zero of 
Jq(x) and the (n + l) s * zero of J^(x). This is graphically demonstrated in figure 15. 
These factors can be obtained conveniently to any degree of accuracy desired from tables 
or books on Bessels function. Averaging these two limits will allow the determination of 


a value of v to be used as a starting point. The companion coordinate /< m can then 
be obtained with the use of equation (27) and the starting value of v just mentioned. 
The modified Newton- Raphson method just described ’’homes in” on the solution ex- 
tremely fast, so fast that an answer can many times be obtained after just three or four 
tries. 
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APPENDIX B 


CALCULATION PROCEDURE 

The procedure used to calculate the U-factors was to initially select a plasma den- 

3 

sity, magnetic field, coil length, and wavelength. A magnetic field of 4. 2x10 gauss was 
selected to correspond with current Lewis Research Center experiments (ref. 1). 
Dimensional frequency was selected to be consistent with (1) the coupling resonant ex- 
pression of reference 2, (2) the previously indicated magnetic field, and (3) the initial 
designated parameters. With this dimensional frequency fixed, a variation of the pair 
(£2, Bq) can be determined. The natural oscillations of a plasma of a given radius R are 
then determined for each (£2, Bq) combination from the simultaneous solution of equa- 
tions (27) and (28). Before the U-factors can be calculated, the radial wave number of 
the forced wave must be found from equations (7) and (8). Hence, for a given coil length, 
a value of v is found for each (£2, Bq) pair. Enough information has been assembled at 
this point to enable the calculation of the U-factors from equations (21) and (36). 
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